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I. INTRODUCTION 


This thesis investigates the use of filtering techniques 
as an electronic counter countermeasure (ECCM) to overcome the 
effects of the optimum jammer operating against coherent phase 
shift key (PSK) and frequency shift key (FSK) receivers. These 
receivers are well-known solutions to the problem of optimum 
detection of known signals in the presence of additive white 
Gaussian noise (WGN) [Ref. 1]. The structure of such receivers 
is shown in Figure A.l. These two receivers represent the 
optimum signal processing algorithms for the discrimination 


of the digital signals S} (t) and sg (t) representing digital 


one" and "zero" logical states respectively. To reduce the 
effectiveness of these optimum receivers, Ref. 2 and Ref. 3 
developed optimum jammers to be used against the PSK and FSK 
coherent receivers. From the resultant jammer waveforms derived 
in these references, this thesis investigates the use of front- 
end filters inserted in these receivers as shown in Figure A.2 
as a method for reducing jammer effectiveness. 

From Ref. 2, the optimum power constrained jammer for 
either the coherent PSK or FSK receiver is a signal propor- 
tional to the correlator signal, shown in Figure A.2 as sq(t). 
This correlator signal (see Ref. 1) is proportional to the 


difference of s} (t) and sg(t), the digital "one" and "zero" 


signals. Though the correlator signal can take different forms 


10 


in these receivers depending on the modulation used, the theory 

used to develop these receivers shows that a jammer waveform, 

proportional to the difference of s (t) and Sp (t), is optimal. 
The jammer, therefore, takes the form of a continuous wave 


(cw) or tone centered at the frequencies of S,(t) and s,(t). 


0 
In the case of the PSK coherent receiver, the jammer is a single 
tone centered at the signaling frequency. For the FSK case, 
the jammer waveform consists of the difference of two tones 


at the frequencies of s(t) and s,(t). The use of tones or 


0 
Cw jammers in angle modulated analog or digital receivers is 
discussed in Ref. 4. 

Given these forms of jammer signals, an obvious choice for 
the front-end filter is a band reject or bandstop filter. The 
filters chosen for analysis in this thesis are the ideal band- 
stop filter and a single zero second order bandstop filter. 
Though the ideal filter is not realizable, the analysis involv- 
ing its use will provide both insight into the problem and 
approximate results before investigating the use of the more 
complex second order filter. For the PSK coherent receiver, 
the bandstop regions of both the ideal and second order filters 
are centered at the frequency of S (t) which is also equal to 
the frequency of So (t). Since in this thesis a filter with 
only a single bandstop region will be analyzed, the bandstop 
region of the ideal and second order filters used in the FSK 
coherent receiver will be centered midway between the frequen- 


cies of S,(t) ana s, (t). 


0 


IL 


The analysis of the two forms will begin with the statis- 
tics of G, the output of the coherent receiver. From Ref. 2, 
G is a conditional Gaussian random variable with conditional 
Statistics depending on whether S4 (t) or Sg (t) was transmitted. 
The analysis will then focus on determining the statistics of 
the random variable G', the output of the receiver with a front- 
end filter as shown in Figure A.2. The mean and variance, 
conditioned on whether s} (t) or Sg (t) was transmitted will 
then be determined using frequency domain analysis techniques. 
Finally, the general analysis will conclude with the determina- 
tion of an expression for the probability of error of the coherent 
digital receiver using a front-end filter to counter the opti- 
mum jammer. Next, the specific frequency domain forms of the 
PSK and FSK signals and jammers will be developed and used to 
calculate the probability of error performance for these re- 
ceivers with and without front-end filters. 

To determine the Fourier transforms of the jammers, this 
thesis will first assume that the jammer waveform is a time 
truncated function. Since in References l, 2, and 3, the 


correlator signal s,(t) is defined over the interval (0,T), 


d 
and the optimum jammer (from Ref. 3) is proportional to Salt); 
the optimum jammer will be treated as defined in the interval 
(0,T) also. Since from a practical standpoint, the jammer is 
a cw signal, the analysis will be modified in the sequel to 
account for this interpretation of the jammer. 


In Chapter II, a development of the mathematical expressions 


in the form of frequency domain equations for the statistics 


T2 


of G' (the output of the coherent receiver of Figure A.2) 

is undertaken. From these statistics, the expression for 

the receiver probability of error is developed. Chapter III 
demonstrates the application of these general expressions to 
the binary PSK (BPSK) coherent receiver problem with both 
ideal and second order front-end filters. Chapter IV presents 
similar analysis to that in Chapter III, applied now to the 
binary FSK (BFSK) coherent receiver. Chapter V presents 
quantitative and graphical analysis of the four receiver/ 
filter combinations, under the assumptions of both the cw and 


time truncated jammer forms. 
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II. ANALYSIS OF DIGITAL COHERENT (RECEIVERS 
WITH FRONT END FILTERS 

In order to investigate the effects of filtering at the 
front-end of digital coherent receivers, it is necessary to 
review the operation and performance of digital coherent 
receivers in the presence of noise and jamming without the 
use of filters to notch out the jammer. As previously pre- 
sented in Ref. 2, the structure of the receiver to be analyzed 


is shown in Figure A.l. In the absence of a jamming signal 


this receiver is optimum for determining whether a "one 
signal or a "zero" signal was transmitted in a given interval 
(0,T), with minimum probability of error [Ref. 1]. 

Following the analysis of Ref. 3, which assumes that the 
jammer is a deterministic waveform unknown to the receiver, 


the input signal r(t) to the receiver front-end is mathematically 


modeled as 


r(t) 2 s,(t) + n(t) + n,(t) 0<t<T, i=0,1 (2.1) 


where the s;(t) (i = 0,1) are used to transmit the digital "one 
and "zero" data, n(t) is a sample function of a white Gaussian 
noise process with power spectral density level of No/2 watts/ 
hertz, and as stated previously no is the jammer waveform. 


From the analysis of Ref. 2, the output of the receiver is 


given by 


14 





I 2 2 
CMM E 3lllsgll - sili ] (2.72) 


where 


T 
(x,y) = f x(t)y(t)at and ||x||^ » (x,x) 
a 


Since G is a (conditional) Gaussian random variable, the sta- 
tistics are completely determined when the conditional mean 


and variance are found. The conditional mean is given by 


E [G| s; transmitted] = (s,,8,) + (n, S4) 
l 2 2 
- Iusti? - Ls, 117: (2.3) 
ve 00, 1 


and the conditional variance is given by 
No 2 
VAR[G|s, (t) transmitted] eae Isa nuuc I (2.4) 


A. FILTERING EFFECTS ON THE MEAN AND VARIANCE 

The effect of placing a filter at the front-end of the 
receiver is now analyzed by evaluating the conditional mean 
and variance of the receiver output. Figure A.2 shows the 
placement of the filter at the receiver front-end. It can be 


seen that G becomes G' where 


15 


1 2 2 
6" = (rs + Hlls [|2-I|ls, | 


(2:59) 


00 


ru = rooms f h(t-a)r(a)da 


and h(t) is the front-end filter impulse response, Since 


r'(t) = s! 


i(t) + n'(t) + Bos i = ORA (2.6) 


it is simple to see that 


ElG'|s, (t) transmitted] = (si,s4) + (nis) 
l 2 2 
- iulsyli?^-Ils li?) (2.7 
1 = 0,1 
where 
S; (t) = S; (t) * n t) 
AE) = “nat E) AE) 
J J 


It is much more convenient to specify filters in the fre- 
quency domain because convolutions in the time domain become 
multiplications in the frequency domain. Thus looking 


specifically at the term (Si Sa? where 


T | 
(sis) = f So SE MI (2.8) 


l6 


Or 


(si,sQ) = E si(t)sa, (t) dt NOE (2.9) 
where 
s = s(t) * pos) 
with 
OR Cas 
pu) = | 
0, elsewhere 


Emsreplacing Seg E by its inverse Fourier transform equiva- 


lent Lu <=> Sa 


(w)), namel 
p y 


Saplt) = LH J 5, (meta (2.10) 


dp = dp 


yields 


' = Ñ ' il i guit 
(ss = poe | S gp wie" “dul dt (2.11) 


Interchanging the order of integration yields 


T 


E T ETE US f 
(e sr) | sape) l J see ad i = OD (20 


= 


and recognizing the fact that 


co 


(f sjte"tat - sic)  i-041 
results in 
(sad e SIN (Daley 
i? ed 2T Ls dp i 5 ` 
Furthermore, since 
! = * 1 = 
s(t) = s; (t) Ip (Ge) jt (os 
Si(w) is given by 
S; (w) = S; (w)H (w) i. = 0,1 
and therefore 
; AA » - : D 
(SiS) Tm f Sgp tw) S;( w) H (-w) dw i = 0,1 (220040) 


“200 
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Replacing sq) with B as inwEguatuoni2c9 7 "theBterrn 


) becomes 


Erg 


00 


(ni,sQ4) = f ni (t) Sq, (t) dt (2.15) 


OO 


and using Equation 2.10 in Equation 2.15 yields 


oo co 


' SS ' JE Jwt 
(ni,sg) = ni (t) l | € p(s)e dw]dt OPEN 


As done n Equations 2.11-2.13, interchanging the order of 


integration and using 


co 


m) = f w at Ro 14) 


— SO 


with 


N' N. 
Y j (w) H (w) 


finally yields 


p= df Sgp (©) Nj (~w) H (-w) dw (2.18) 
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The conditional mean of g', given that s; (t) is transmitted, 


is 


; l 
l = — 
E [G s; (t) transmitted] F | S 


=O 


1 
tae J Se eh LS 


«= 00 


2 
| 


l 2 
* iis LH? - Ls, (2.19) 


i = 0,1 


The variance of g' given that s(t) is transmitted can be 


shown to take on the form 
VAR[G'|s; transmitted] = A (2.20) 
The right hand side of this equation becomes 


Eia O 


00 00 


= a SOS _J aa uM (2 2) 
Equation 2.21 can now be written as 
E[(n'(t),s,(t)) 7] 
Sa 
a E E[n' (t) n' (t) sap (t) sap (7) dtar (2192109 
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Recognizing that the expected value operation in the integral 
is the autocorrelation function of the filtered white Gaussian 


noise, Ri. (t-t), Equation 2.21 can be written as 


= f a) Rat (£71) $4, (t) sq (t) dtdt (2.23) 


Again, it is much more convenient to continue the analysis 
in the frequency domain. Using the inverse Fourier transform 
for the frequency domain representation of R, 1 (t-7) and Sub: 


stituting into Equation 2.23 yields 


1 ju (£-1) 
- [f JG f Sq 001677" àdlsa (E) s, (1) dtár 


(2.24) 


where a! is defined by Equation 2.9. Interchanging the 


order of integration, Equation 2.24 becomes 


2 


El (n',sg) ] 


o oo oo 


l 
= l 3, 1 (w) du f E 


It is easily recognized that 
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"e jwt 
f sap (67^ dt - Saj(-e) = Sã (o) 


— (O) 


and 


ù = JOT = 
f Sap 7$ dt Sap tt) 


00 


Therefore in the frequency domain, VAR(G'|s; transmitted) 


becomes 


eo 


ll 2 
A f Sp tu) IS gp tu) | dw (2.26) 


== 00 


VARIG '|s, (t) transmitted] = 


Finally since n'(t) = n(t) * h(t), it is well known that 
2 
SUUM S (w) |H (w) | 


Since n(t) is a sample function of a random process that was 


assumed to have power spectral density S (w) = N9/2, for all w, 


N co 
' , 0 2 2 
Var [G |s; (t) transmitted] = = E |H (w) | IS ¿y (0) | dw (29285 


where acEME I (2) [418 4,00 |^ is an even symmetric 
function has been used to reduce che mes Pol intao om 


to half the real line. 
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In summary, G' is a (conditional) Gaussian random variable 


with conditional mean and variance given by 


E[G'|s; (t) transmitted] 


co 


2 = * Sap (9) $4 (70) B (C70) do 


co 


il 1 2 2 
ot TA f oe ei +5[115,| | -{{s, || ] 


— OO 


IPSO] (2.28) 


and 


VAR[G'|s; (t) transmitted] 


co 


J l|*]sq ()|^àe — i - 0,2 (2.29) 


o 
2T 
B. PROBABILITY OF ERROR PERFORMANCE 

The analysis thus far has been concerned with determining 
the statistics of the random variable G', the output of the 
coherent receiver. The next step is to determine the proba- 
bility of error from the decision process, where the amplitude 
of G' is compared with a threshold y. Referring to Figure A.2, 
if the value of G' is greater than y, the receiver decides that 
a "one" or s, (t) was transmitted, and conversely if the value 
of G' is less than y the receiver decides that a "zero" or 
Sg (t) was transmitted. The quantity, y, is determined in the 


derivation of the optimum receiver for detecting known signals 
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in additive white noise [see Ref. 1]. That analysis determines 


the threshold of the correlator (coherent) receiver as 


N T 
Y 7" —Ginne*z| ya (ea eee (2.30) 
0 


However the receiver of Figure A.2 treats the second term of 
Equation 2.30 as a bias, thus the threshold becomes 
No 

A cR ne (2310) 
where Ng/2 is the power spectral density level of the white 
Gaussian noise process. Here, n is determined by the specific 
decision rule applied. In this case, the receiver strategy 
is to detect s, (t) Or sg (t) with minimum probability of error, 


so from Ref. 1 


P{s. (t) was transmitted} 
"m B (2.990 
Pisg(t) was transmitted] 


All further analysis in this thesis will assume that the 
probability that s; (t) was transmitted is equal to that of 


sg (t) being transmitted, so the decision threshold is 
g&n 1 = 0 (2 75959 


Now, the receiver can make two types of errors.  Specif- 


ically, deciding that So tt) was transmitted when in fact S, (t) 
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was transmitted (defined as error 1), and deciding that s, (t) 


1“ 
was transmitted when So (t) was actually transmitted (defined 

as error 2). The receiver decides that So (t?) was transmitted 
when G' is less than the threshold y, and that sy (t) was trans- 
E ted when G' is greater than y. So, the probability of error 
l occurring is the probability that G' is less than y when s} (t) 
is transmitted. Since G' is a Gaussian random variable with 
conditional probability desnity functions dependent on whether 
sy (t) or sg (€) was transmitted 

(g'-m,)* 


2N 
e J ag! (234) 





Perror 1 


where m, - E[G'|s, (t) was transmitted] and Me VAR[G' | s, (t) 
was transmitted]. By letting x = gm /VV 1 Equation 2.34 
becomes the integral of the standard Gaussian probability 
density function also known as the error function (erf(x)) 


of the form 


P dx (2.35) 
error 1I 


—00 /27 


The probability of error 2 occurring is the probability that 





G' is greater than y when sg (t) is transmitted. Therefore 
-(g'-m P 
P = i E E : J dg' (2.36) 
error 2 E 
Y Penna 
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where 


My = E[G'|s, (t) transmitted 


Using a similar change of variables as in Equation 2.35 


yields 


00 


il -x^/2 
Perror Ql p x 
-mg/ uS y 2T 


dx (22370 


where y is zero from Equation 2.33. Finally, the average 
probability of error becomes 
(250959 


P = pP + (l-p) P 


e error l error 2 


where p = Pis, (t) was transmitted) and 1-p = Piso (t) was 
transmitted). Since p = 1/2 and substituting Equations 2.35 


and 2.37 into 2.38, the average probability of error is 


“mo LG 
my, / ps 








2 
Po = zi í Lone He a 
< —0 /27 
ELE! -x?/2 
4 f e dx] (2:39) 
-mo/ "v, /27 


TITR PHASE SHIFT KEY COHERENT RECEIVER 
WITH A FRONT-END FILTER 

A. PHASE SHIFT KEY WITH AN IDEAL FRONT-END FILTER 

Using the equations for the mean and variance of the con- 
ditional Gaussian random variable G' (Equations 2.28 and 2.29), 
along with the probability of error equation (Equation 2.39), 
it is possible to investigate the effect of a front-end filter 
on the probability of error performance of the coherent receiver. 
The receiver analyzed in this case is the optimum receiver for 
binary phase shift key modulated signals with an ideal band- 
stop filter in the front-end. The filter will be characterized 
as having a gain of unity over all frequencies, with the excep- 


tion of the regions from w = 209 -B/2 to w = tw, +B/2 where 


0 
the gain is zero. The bandwidth of the filter is B radians/ 


second and the center frequency is uo, radians/second. Although 


0 
this filter is not realizable, it serves both as a good approxi- 
mation to a higher order bandstop filter, and as a simple 
example so as to gain insight to the receiver performance. 


Using Figure A.2 as a reference, the generic coherent 


digital receiver can be made into a BPSK receiver by setting 
sq) = 2A cos Wat OS I C3215) 


and modeling the input PSK signals as 
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Ss, (t) = A cos (wt + 9; (t)) O <S E <p io (35:2) 


where 
T— for ^1 - © 
9 
rorem 
Thus s. (t) (i = 0,1) can be expressed as 





sy (€) = -A cos Wgt 
ee ( 32899) 

s (t) = A cos Wot. 

The jammer waveform is set to 
n; (t) = E ¥ 2/72 cos Wet Quo c m (3.4) 


since in Ref. 2 it was shown to be optimum against BPSK trans- 
mission. Here E ls the jammer output power. Finally the 
noise input into the receiver of Figure A.2 is a white 
Gaussian random process with power spectral density level of 
Ny/2 watts/hertz. 
l. Calculation of the Variance 
From Equation 2.27, the variance of G' conditioned on 


s; (t) being transmitted is given by 
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VAR[G' | s, (t) transmitted] 


00 


N 
= = J I (a) |^] Sa. (2) | "au i=0,1 (3.5) 


Since the front-end filter is ideal, Equation 3.5 can be 


modified to become 


VAR[G'|s, (t) transmitted] 


- O, f js Ola - J | s to | au (376) 
2T 9 dp i -B/2 dp 


since owe is unity except over the bandstop region. 


2 : 
In order to calculate Sap tu) | , since sap (t) = s,(t) p(t) 
where 
ILL ST 
p(t) = 
O elsewhere 
then 
Sap E) = 2A p(t) cos ppt (30999) 


The Fourier transform identity f(t)-g(t) <=> (1/27) xF(w) * G(w) 


will be used in order to calculate Sap 9 - Since P(w) is 
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Payment! Saar/2) ene E (3.8) 
and 
Sq (w) = 2An [6 (709) + $ (009) ] 


Sapo? is the convolution of these Fourier transforms, result- 


ing in 
T 7j (»-09) 1/2 
Sap tt) = AT [Sa (n709) 5€ 
T -j (w+w g) T/2 
+ Sa (wtw,) xe ] (3:298 
0 
From this, IS45 0 |* is obtained by multiplying Spe by 
Sg (-w) which results in 
P 

2 2 E T E T 
Sap (1)! = A'T [Saí w ugly Sa (wtu) (3.10) 
where the product of Sa (w+w,) T/2 and Sa (0709) T/2 has been 
ignored because it is essentially zero for practical values 
Of wo anden: 

Since the variance consists of an integral fromw= 0 


to infinity, for large Wo the contribution of the Sa (u+w,) T/2 
term is negligible. Therefore, substituting into Equation 3.6 


yields 
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Var [G'| s; (t) transmitted] 


N o9 *B/2 
2,12 o 
rer 9 TES da f salo-wp)5dwe] ^ (3.11) 
2T 0 0'2 02 
wg 7B/2 


i = 0,1 


In this form Equation 2.11 is unmanageable due to its complexity 
and the number of variables present. To simplify, the follow- 


ing substitutions will be made: 
l. x= (0709) T/2, 2/Tdx = dw 


2. for wu = 0, x = “Wy T/2 which for large values of WoT 


becomes approximately -o 
S. for = o, X = 0 
4. for wu = Wo no (2 Xx = +BT/4 


Equation 3.11 now becomes 


VAR[G'|s; (t) transmitted] 





Aris o o BT/4 
x Lf Seme acs = f Sa (x) dx] (02) 
T^ -e -BT/4 
i= Ok 
| 2 
Finally the Sa(x) is a symmetric function about x = O and 


f dd = 1/2 from math tables. Thus the variance of C' 
0 | 


3d 


can be written as 
VAR[G'[s; (t) transmitted] 


2 BT/4 


= PN AR Za - if Sa^ (x) dx] Al! (So 
0 e 0 


2. Calculation of the Conditional Mean 
The conditional mean of G' is given by Equation 2.7 


and shown here for convenience, 


E[G'|s; (t) transmitted] 


» í 2 l DM 2 M 

= iss (ni, S4) HEURE lls) 1^3 i= o 
and consists of three terms. The frequency domain forms of 
(Si 7S) and ea are given by Equations 2.14 and 2.18 
respectively. The term 50] 15) 115 -| Is, 1173 (recall that 


T 
LESS = di Su represents the difference between the 
0 
energy contained in sg (t) and s, (t). Since our signals are 
equal energy signals, the term in question becomes zero. Thus 
the conditional mean becomes 
Jd -is Om (3.14) 


E[G'|s; (t) transmitted] -  (si,s4) + (ni,s 


a? as! 


a? and recalling that 


H(w) is an ideal filter, Equation 2.18 can be modified to 


Looking first at the term a 


52 


' = 0 = 
E zu J Sao N w) dw 
wgtB/2 -wg +B/2 
= Al S.  (0)N. (-o0) du - S 7 _(0)N, (-w)dw] 
0 =B/2 dp J -w -B/2 dp j 
0 
(SS) 
in order to evaluate Equation 3.15, the product ERAS (0) 


must be calculated using nj C6) as given by Equation 3.4. 
However, since the model of the coherent receiver in Figure 
A.2 contains a finite time integrator, put is effectively 
truncated to the interval 0 < t < T, or equivalently, the 


JP 
p(t) is as previously defined. Thus, Se Un. becomes 


new term n. (t) replaces ata where Mpio = DS So and 


H — T 
N. (0) = E voy T > 


T -j (w-wg)} T/2 -j (tog) T/2 
[Sa (0709) 5e t Sa (w+w,) T/2e ] (3-16) 


Mültiplying me and Sap e? yields four terms: two Sa terms 


centered at Wo and “Wo 


Or Sa (otog) T/2 and sa (0-09) T/2. It can easily be seen that 


and two terms that involve the product 


for large w, this product is essentially zero. Substituting 


0 
into Equation 3.15 it is also easy to see that the contribution 
of the Sa(wtw,) T/2 term in the integral with limits of +wy+B/2 


is negligible, so that 8 Ba) becomes 
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j d ¿TZ nj 
= wgtT/2 
1 sa ) Wy = J Sa^ (ut ) Lau) 
0" 2 0-2 
= 00 w9-B/2 
m “Wy +B/2 
2 T 2 T 
+ ( E Sa (wtw) 57 f Sa (w+w¿) 7dw) ] 


- 7-0g-B/2 


(328 


Recognizing that the integrals of Sd(wtw,)T/2 and SA(w-w)) 7/2 


are equal, performing the same change of variables as in 


Equation 3.11, mE becomes 


= MD 271 
(ni, Sq) = = Paj ve a 
E BT/4 > 
[ f sa (x)dx - Í Sa (x)dx] 
— o0 -BT/4 
Finally, recognizing symmetry about x - 0 and usihg 


f saíx)dx = 1/2, we obtain 


BT /4 


: E 77s 2,70 _ 2 
(n5, Sq) = AT Pp: /2/T = ley E Sa“ (x) dx] 


(3. 180 


(3 19) 


Looking now at the term (Si rS)? recalling again that 


H(w) is an ideal filter, (s; SQ) can be written as 
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ogtB/2 


oo 


ros = =| f Sgp (0) 8, (79) do - 


S. (w)S. (w) dw 
= 4. 7B/2 E EE 


709g-B/2 


- it S, (w)S, (-w) dw] (3.20) 
-ug-B/2 EDS 


q 


The signal models for sg (t) and s, (t) given by Equation 3.3 


can be Fourier transformed to yield 


-j (w- 2 
| AT NR IDA 
So (w) = - 5- ISa(o709) 5€ 
T -j (tog) T/2 
+ Sa(wtw,) ze ] (27215 
0. 2 
and 
nm T Z 
_ AT E 0 
91 (w) = > 18 a(w wg?) € 
-j (w+wg) T/2 
+ Salwtw,) xe (32 2) 
055 32 
To calculate the specific case of (Sj sS)? iu St EGET 


the product of CE This yields four terms as in 


Re where two are Sa terms centered at two and two are 


products of Sa (w+w,) T/2 and Sa (w-w5)T/2. The latter products 


can be neglected for large w Also, with respect to Equation 


0" 
ENDO the contribution of the Sa(w+w,) T/2 term to the integral 


whose limits are tog *B/2 is negligible. Therefore, substituting 


Sap 2) 9g (79) into Equation 3.20 yields 
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222 oo 





(sos) = 32 m Sa (w-wp)T/2 du 
og *B/? = 

S: n sa? (u-wp) du) + ( EN ) du 

w.=B/2 2 0 

0 
-w -B/2 
9 2 T 

- f Sa” (w+w,) 5 du)] (3 28) 

-09-B/? 


2 
It is easily seen that the integrals involving Sa (utw,) T/2 
are equal to the integrals involving the Sé (w-wp) T/2 term. 
Making use of this observation and the same change of variables 


used in Equation 3.9, it is easily seen that 


' AA 2 à ii BT/4 
(59:94). 7 -A T «|. f saix)ax - f | saá G0 dx] (3 4265 
-o -BT/4 


2 
and by using the symmetry of the Sa (x) function, (sor Sa! 


becomes 


pur 
(S015) Gee EI d Sa^ (x) dx] (328) 


By similar analysis it can be shown that 


BT/4 
(st,s,) = a'r Z(t - f sa? (x)dx] (3.26) 
0 
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In general terms, (s;,Ss¿) (1 = 0,1) can be written as 


| BT/4 
Pee T | GS we) i-0,1 om 


so that the conditional mean of G' can be expressed as 


E[G'|s; (t) transmitted] 


BT/4 


= (ATVB, V27T ee ) aer) (SCD - fsa? (x) ax) } (3.28) 
0 


1= 0,1 


See RrODability Of Error Calculation 
With the conditional mean and variance of G' given by 
Equations 3.28 and 3.13, it is now possible to compute the 
probability of error performance of the BPSK coherent receiver 
with an ideal front-end filter. Assuming the probability of 
the transmitter sending a "one" or a "zero" signal is the same 
(P(si (t)) = P(sg(t)) = 1/2)), Equation 2.39 shows that the 


average probability of error Po is 


00 2 m/v, 2 
[ f e t f e * /*ax] (E 


-mo/ "v. /27 =% 


ho 
II 
N| FP 


where the threshold setting given by Equation 2.33 has 


been utilized. 
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In order to compute Por all that is needed is the 
specific conditional mean divided by the square root of the 
conditional variance and using this ratio in the limit of 
the appropriate error function integral. First, for a "zero" 


transmitted, 


2 
x RE X (3.30) 
BT/4 
[5 - f sa? (x) dx) 
0 


BT/4 > 
[AT/P 5/2/€ -A 2m] [É(E - f Sa“ (x) dx) | 
0 


| 

als 

Q e 
NO 


al 


A TN, 


2 
Recognizing that A T/2 = E, the average energy per bit and 


substituting this in Equation 3.30 one obtains 








> 7 BT/4 , 
Em [/4gPh; - 2bIU- C LM ad (SDN 
0 J E 0 
a (39598 
V BT/4 
2 2 
3 2EN, f (5 - a Sa^ (x) dx) 
Finally 
BT/4 
0 A HÀ 2 
p /2E/NQ(12-YP,; 7E) T(5 - |  sa® (x) dx (3.32) 
V 0 
g 
and similarly 
-my BT/4 > 
sae = AYN Le ;/E) 3 - f Sa (x) dx) (3 08) 
d. 0 
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By using Equations 3.32 and 3.33 as the limits of integration 
ataron 3729, the probability of error performance can now 
be determined versus E/Ng, the signal to noise ratio; Paj? 
the jammer to signal ratio; and BT, the relationship of the 
filter bandwidth to the inverse of the signal bandwidth. 

Several items of note can be observed in the calcula- 
tion of P_. First, as BT approaches zero, the term 

A | | 

2/7"T(m/2 - f Sa (x)dx) becomes unity and Equation 3.29 
becomes e receiver probability of error under optimum jamming 
[See Ref. 2] for BPSK modulation. Furthermore, when both BT 
and Eum the jammer to signal ratio, go to zero, Equation 
3.29 becomes the well-known result for performance of a 


coherent BPSK receiver operating in additive white Gaussian 


noise [see Ref. 1]. 


B. COHERENT PSK USING A SECOND ORDER FRONT-END FILTER 
Instead of choosing a filter transfer function that 
approaches the ideal case, the effect of a second order front- 
end filter on the probability of error performance of the 
coherent BPSK receiver will be investigated. The general 
coherent receiver structure becomes the BPSK coherent receiver 
by using Equations 3.1 and 3.3 in place of sq) and s; (t) 
(i = 0,1) respectively. The filter will be characterized as 
having a single notch at the center frequency of the jammer. 
The filter transfer function H(s), is derived from a first 
order, low pass prototype given by 


(3.34) 


H (s) = == 
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To transform this into a bandstop filter centered at wg’ the 
frequency transformation 
S (o ED) 
a —_ (3.35) 
S FW Wr 


will be substituted in Equation 3.34. Here, Wry and o, repre- 


sent the upper and lower 3 dB frequencies of the bandstop 


region. Before substitution, it can be seen that for wr 


approximately equal to w or equivalently w, -w. small com- 


U’ U L 

; 2 
yt’ Vyr and Wi, TW, are approximately equal to Wo 
and B respectively. Note that B is the half power bandwidth 


pared eon. 


of the bandstop region. Thus, Equation 3.35 becomes 


s mU (3.36) 


s* +06 


Substituting finally into Equation 3.34 and also setting s - jo 
in order to arrive at the frequency domain form, the filter 


transfer function is 


H(w) = = aa (IIS 
-w  juB 
This filter will prove somewhat cumbersome to work with when 
put into Equations 2.147 2.18 ane 22H uinererore almere 
convenient form of this filter will be sought later. 

Looking at Equations 27140192 0155900012 EN SAUS Es 


2 y 
Sap (9) $1 C79) ; Sap (9), CCo) and IS gp (0) | (calculated in the 
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dersvatron of Equations 3.10, 3.18, 3.23) are Foureir trans- 
forms of a time limited signal multiplied by a cos (wot). The 


time domain forms of these products are of the general form 


c(t) = a(t) cos (wt) UNE. ST (95 58) 
which is essentially an "AM" modulated signal. This is not 
surprising because each signal in the coherent BPSK receiver 
model is of constant phase over a bit interval (0 « t < T). 
Reference 5 shows that the effect of filtering a modulated 
cosine wave can be determined by calculating the effect of the 
modulation being passed through an equivalent baseband filter. 
The baseband equivalent of the bandstop filter will be in this 
case a high pass filter. 

The equivalent high pass filter model can be derived by 


applyinc the transformation 


Q 
E, M (3.39) 
S 
to the low pass prototype of Equation 3.34, and letting i 
the upper cutoff frequency be equal to B/2. The resulting 
filter in the frequency domain is 


B/2 + Jw 


When Equation 3.37 is evaluated under the assumptions of 


large Wo and small B compared to o: the filter of Equation 
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3.40 can be seen as an equivalent, frequency shifted version 
of the filter of Equation 3.37. Using the filter of Equation 
3.40 and just the "modulation" spectra of the products, the 
effect of a second order bandstop filter in the front-end of 
the BPSK receiver can be evaluated with respect to the proba- 
bility of error performance. 
l. Calculation of the Conditional Mean 

The conditional mean of G' is given by Equation 2.7 
where n4 (Sq) and (Si Sa)? are given by Equations 2.14 and 
2.18, with the added fact that ai 118,115 is zero. 
There is now enough information to determine the conditional 


mean of G' namely 
' : e t t > = 
E [G EN transmit teda (n;,sg) + (si;Sq) OA 


Taking each term individually, ON will be calculated 


first. Equation 2.18, repeated here for clarity, 


1 o 
f NI a 


ae OD 


) ES 


8 


shows that the product Sdp ON and H(-w) will be needed. 


The product S ey was previously found (in the derivation 


d 
ot. Equation 3 Es be 


2 
o NE: T 
S ap (#) N. Co) = VP aj” / TISE (o wg) 7 +tSalwtwg) z] (3.41) 
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Before multiplying this by H(-w) and substituting into Equa- 
tion 2.18, rt can be seen that Equation 3.41 can be written 


in the form 


a NE (0) 


= ATVP j “27T Isa” (w) 5 * m[8(u-ug) *6 (atu) ] D 


which in the time domain is of the form of Equation 3.38. 
Thus the baseband equivalents of ELI and H(-w) can be 


substituted into Equation 2.18. Therefore (ni (Sa) becomes 


eo 


(n;,Sg) = IATP jT P A da (9 49) 


where Tsa(wT/2) is the modulation spectrum and H'(-w) is the 
equivalent high pass filter Finally, substituting Equation 


3.40 into Equation 3.45 yields 


o 


=- LarvVP 277 : peer) 
(ajsa) = gr ATV/P../2/7 J Isa” (wT/2)- 35877 dw (3.44) 


Now taking advantage of the simplified expressions for the 
filter, the "modulation" term, and the respective time domain 
forms, Equation 3.44 can be evaluated using Parseval's theorem. 


Parseval's theorem states that if £, (8) and f. (t) are time 


2 


functions with Fourier transforms F4 (0) and F, (w) then 
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00 il oo 
f £,(t)f,(t)at = 3 [ Fi(w)Fj(-u) do 


By letting -jw/(-Jw+B/2) equal Fo (cw), it is easy to show the 
impulse response of the baseband filter h'(t) is 


h'(t) = 6(t) ons 


u(t) (3.45) 
Before using Parseval's theorem, the inverse Fourier transform 
of T(SE(w) T/2 must be found. This is easily determined from 


tables to be 


e3| ct 
+ 
— 
i 
FH 
[^ 
ct 
|^ 
e 


EISE TRUE | (3.46) 
[B 


| 
H3| ct 
© 
|^ 
d 
|^ 
= 


Putting S ) into the equivalent time domain integral form 


d 


using Parseval's theorem, yields 


Bi e) de EA 


T 
; à SE n _t .B 
(ni,sq) = AT /B.. /2/T Jl (a1) (56€) Ze 
due to the fact that the filter is causal (only defined for 
t > 0) and that r liTSaĝĵw)T/2) is defined from -T < t< T. 
Evaluating Equation 3.47 consists simply of the task of evalu- 
ating the sum of four integrals with limits from zero to T. 


It is easily demonstrated that 
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uH — 2 — -BT/2 
DES S AT Ps /2/T [sm e ) ] (3.48) 


a) 


The second term that must be evaluated is (Si;SQq); 


(i = 0,1) in order to be able to determine the conditional mean 
of G'. From Equation 2.14 repeated here for convenience, 
1 00 
' = Mu m = | = 
| Sap (09, l w)H(-w)dw i = 0,1 


«=> 00 


It is apparent that the products of Sm and Si (70) (i = 0,1) 
have been determined in the derivation of Equation 3.24 and 
can easily be put in the form of Equation 3.42. Therefore, 
the baseband equivalents of the filter and the product are 


substituted into Equation 2.14 which yields 


oo 


DM 
A => 
I f TSa” (WT/2) =30+877 


(S9: S4! = zm du (3.49) 


It is easily seen that the integral in this equation is of 
the same form as Equation 3.44, thus Parseval's theorem can 


be applied directly and for the specific case of (SQ (Sa) yields 


T 
T ou -t _B_-BT/2 
(Sg) = “AT "Nc RR E S t eee (550) 


Here the limits of the integral are 0 and T because the filter 


impulse response h'(t) is causal. Since the integral has 


45 


already been evaluated, we obtain 


(Suo -A^T[; e (1 ce PI/2] ENS 


and similarly 


= (3.52) 


, A 2 2 B 
(s ; Sq) = A Tig e 
Finally, the conditional mean of G' under the assump- 


tion that a "one" or a "zero" signal was transmitted is 


E[G'|s, transmitted] 


= [aT VP . V27T 4 -1 2^5 nea eee (3:358) 
nj BT 
JUNI 
2. Calculation of the Conditional Varwance 


In order to calculate the conditional variance of G', 
it will be more convenient to use Equation 2.27 due to the fact 
that Parseval's theorem can be applied to its evaluation. 


The conditional variance of G' is 


VAR[G'|s, (t) transmitted] 


Z 


Co 


= p f Ew) [Sgp |f dw i= 07 (3.54) 


== 00 


The term |S o tu |? has been calculated and is given by Equation 


d 


3.10, and as shown here, it is of a form similar to Equation 
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3.41, namely 


2 
[S4,,(0) |^ - A°T* [satu-w9)T/2+ sa” (wtwg) 7/21 (3.55) 
Therefore the baseband equivalent of [Sgy tw) 1%, that is 
2 
São (w) |* = 2AºTº [Sa(wT/2)] GE 


will be used along with the baseband equivalent filter H' (w). 
Eul’ is calculated by multiplying H' (w) by H' (-w), producing 
the following result 

2 


JH’ (w) |f = NEVER 57) 
(B/2)" +u 


Thus the conditional variance of G' becomes 


N 2A 7 
VAR[G'|s. transmitted] = EP V 
aL 4T 
o 2 
[ T$a^(wT/2) [——55——,]ds i - 0,1 (3.58) 
= (B72) Fu 


| cuoscotucing Equations 3.56 and 3.57 into 3.54. 

Use of Parseval's theorem in order to compute the 
conditional variance of G' requires the inverse Fourier trans- 
2 2 PELIS ! ! 

forms of TSa(uT/2) and o /(H2)" +w°). The first inverse 


transform has already been determined and is 
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Po rs IE (3.59) 


— 
I 
F3| rt 
at 
a 
© 
|^ 
ES 
RA 
H 


Next, the inverse Fourier transform of Im | must be deter- 
mined. Since lH* (w) |? was obtained by multiplying H'(w) by 


H'(-w), Pru a TO can be written in the following form 
Flat = Fil + Firm) (3.60) 


Now, using the Fourier transform identity EH > hos 
Equation 3.60 becomes 


B -Bt/2 


ppa) | S IDÉES, -3€ u(t)] * [8 (-e) Be 7 


u(-t)] 


(3.619 


where the necessary inverse Fourier transform has been deter- 
mined in Equation 3.45. Substituting this into the convolution 


integral yields 


co 


Fire |) = de 


— OO) 


BB y (2 ] 


[8 (-x) - PCIA ar (3.62) 


and performing the convolution specified in Equation 3.62, 


200. 
the inverse Fourier transform of |H'(w)|” is 
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EIE LL E -F eBltl/2 SEE 


Using Parseval's theorem, the conditional variance of 


G' becomes 


VAR[G'| s; transmitted] = APTN, 


0 i 
"ED f G41) (8 (te) -ZeP/*yae +f (FE 41) (6 (e) - De PU? 
0 


) at] 


1. 205 T (3.64) 


The evaluation of Equation 3.64 is fairly straightforward, 

with the exception that care must be taken when integrating 

the 6(t) over -T to zero and from zero to T. There is only one 
singularity present at t = 0 with an area of one, so that each 

integral of ó(t) can be treated as being equal to 1/2. There- 


fore, the conditional variance of G' is 


VAR [G' | s; (€) transmitted] 


a 2 2 -BT/2 Nas 
= A IN, lg (t e ) ] i= 0,1 (3-695) 
2 PoPa bilitv of Error Calculation 
With the conditional mean and variance of G' specified 
by Equations 3.53 and 3.64, it is now possible to compute the 
probability of error performance of the BPSK coherent receiver 


with a second order front-end filter. Assuming the probability 


49 


of the transmitter sending a "one" or a "zero" signal are 
equal (P (sy (t)) - P(sg(t)) = 1/2)), Equation 2.39 shows that 


the average probability of error Pa is 


oo 2 
Ea = >| j ei /2ax 
-mg/ MC YAT 
"Ino 
A 
+ f —— e dx] (3.66) 
= 00 y2T 


where the threshold setting given by Equation 2.33 has been 
utilized. 

TO compute si all that is needed is the specific 
conditional mean divided by the square root of the conditional 
variance and using this ratio in the limit of the appropriate 
error function. First, for a zero transmitted 


[A^T - AT/P.- (271) (1 Sen? “ayer 
= nj e (3.67) 


vam (1 -e 9/2) 2/r 


| 
ad 
Q O 


Recognizing that A^T/2 = E, the average energy per bit, and 


substituting this in Equation Brocas 


= 
-— = V2E/N, - / (2E/N9) (PL; /8) ( V (1 -e7B7/2)2/8r) 


g (3260) 


Finally 
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= /2E/N, (1 = a4, a EE E (3:469) 


AÊ 


and similarly 


- - VEEANSUPLTE 1) VQ EI (3.70) 


i 
E 
Q F2 


The probability of error performance can now be determined 
versus E/Ng, the signal to noise ratio; Pa3j4 Er the jammer to 
signal ratio; and BT, the relationship of the filter bandwidth 
to the inverse of the signal bandwidth, by using Equations 3.69 
ume 3.70 as the limits of integration in Equation 3.66. 

To evaluate the probability of error performance of 
the BPSK coherent receiver as the filter bandwidth approaches 
zero, L'Hopital's rule must be utilized. Applying this rule 
eemzaquations 3.53 and 3.65 and then SUBSEA VEN, the results 
into Equation 3.67 shows that, as the filter bandwidth (multi- 
plied by the bit length T) goes to zero, the probability of 
error performance of the filtered coherent BPSK receiver 
becomes that of the coherent BPSK receiver in the presence of 
noise and jamming derived in Ref. 2. Furthermore when the 
jammer to signal ratio goes to zero, the probability of error 
performance is the same as that of a coherent BPSK receiver in 
the presence of noise along as derived in Ref. l. 

As one more check, when the filter bandwidth approaches 


infinity both Equations 3.69 and 3.70 become zero and the 


S 


probability of error is one half. This exhibits similar 
behavior in the limits of BT as in the calculations performed 
using the ideal filter. In a follow-on chapter, the behavior 
of the BPSK receiver with either type of front-end filter 


will be compared over the complete range of variables. 
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IV. FREQUENCY SHIFT KEY COHERENT RECEIVER 
WITH A FRONT-END FILTER 

A. FSK WITH A SINGLE IDEAL FILTER 

As done with the BPSK receiver, the FSK receiver will be 
evaluated using an ideal front-end filter, first to gain in- 
sight into the problem and second to obtain an approximation 
to the performance achieved using a higher order filter. In 
the first case, the front-end filter will be modelled by an 
ideal bandstop filter with a center frequency to radians/ 
second and the bandstop region defined from £o. -B/2 to tw +B/2 
where B is the bandwidth in radians/second. 

In order to obtain a coherent FSK receiver from the 


Ene rico coherent receiver of Pigure A.2, sqQ(t) is set to 


uU. Ae A E 
sq (t) E 2A sinz(w, wa) t sin5(o,tug)t e 1) 


0 SAL 


The input FSK signals are 
sy (t) = Acosw,t OS eae 


and 


sg (t) DL Cos ire Ds 


0 


^ 
ct 
| A 
r3 


(472) 
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The optimum jammer derived in Ref. 2 for use against the 


coherent FSK receiver (without a front-end filter) is 


2 ST 4% ein Em - men 
n, (t) - NE m sin 5 (91 vg)t sins(w two) t (4.3) 


0 MESE 


Finally, the additive noise is modeled as a white Gaussian 
noise process with a power spectral density level of Ng/2 
watts/hertz. 

The first step in determining the probability of error 
performance (as previously done) involves determining the 
statistics of the random variable G' at the output of the 
receiver. Since it can be shown that G' is a conditional 


Gaussian random variable, all that is needed are the conditional 


mean and variance under the assumption that a "one" or a "zero" 
are transmitted. Finally, the probability of error performance 
will be calculated by applying these statistics to Equation 
23037 
l. Calculation of the Condrtrons SN 
In order to determine the conditional mean, Equation 


2.7 becomes the starting point. Equation 2.7 is repeated here 


for convenience, that is 
E[G|s.(t) transmitted] = + (Mis ) «itis n -||s m 
al i d js 2 0 JL 


It has already been shown in the PSK system analysis that 


for equal duration, equal amplitude signals, the term 
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yeg ls, 115 is zero. The frequency domain equations 


for (Si Sa) and DE are given by Equations 2.14 and 2.18 


a) 
respectively.  Modifying Equation 2.14 in order to account 


for the fact that H(w) is an ideal filter, (S: S4) becomes 


oo 


: o T 
(Sir Sa) = ant J RA MU = 
w +B/2 -w -B/2 
S, (=w) S4 gm REM 
-w -= 7B/2 P -w.-B/2 
des (4.4) 


In order to calculate (Si rS)? Equation 4.4 shows that the 
product of S. (=w) and Seg is needed. Since s;(t) is defined 


over 0 « t « T, it can be shown that S, (w) given by 


PEN 2 
AT A 
S. (w) = > [Sa (w-=w,)>38 
-j (w+w;) T/2 
+ Sa (wtw. ) 5€ ] TS 0t (4.5) 
En e) is determined to be 
-j (w-w,)T/2 -j (wtw,)T/2 
_ AT D i T 1 
Sgp = — [Sa (w 04)35e +Sa (w+w] ) 5€ 
- -j(w-ug) T/2 T -j(u-09) T/2 
- Sa (709) 5€ -Sa (w+w,) 5€ ] (4.6) 
Before investigating the specific case (SQ Sq): ae te 


will be convenient to look at the relationship between the 
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center frequency of the filter and the signal and jammer fre- 
quencies. Since the filter contains a single bandstop region 
and the jammer has a spectrum similar to that of eae the 
filter center frequency should be located between the jammer 
frequencies in order to notch out the jammer. Also, the 
filter should have a bandwidth sufficient to cover both w 
and Q4. Although the midpoint of wo and Wy is given by 


(0g to4)/2, it will prove to be much more convenient to define 


w. and w, in terms of Wo’ that is 


jJ 0 
Q4 = d + wa 
(4.7) 
20 Mica oad 
where Qa 7 (0, 709) /2 and Qo = (0, 09) /2. 


Substituting Equation 4.7 for Wy and Wo and multiplying 


St (O by Sg (C70) yields 


Jum -juqT 


T E E T 
[Sa (w=w tug) 7Sa (o7 t04) 5e +Sa (uto 704) 58a (oto +u .) 5e 


2 T 2 Tr 
-Sa (o70 tu) » -Sa (WW 70 51 (4.8) 


Before the substitution of Equation 4.8 into Equation 
4.4 is made, some simplification can be achieved. The contribu- 


2 
tion of the Sa(x) functions with arguments of (wtw i04), to 
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the integrals with limits of tw tB/2 is negligible for large 
Wo and small B compared to Woe This also holds true for the 
products Sa (wtw +wa)T/2:Sa(wt+w -w,) T/2 and Sa (w-w tw) T/2 

o (070 704) T/2. Therefore substituting Equation 4.8 into 


Equation 4.4 yields 


Ds Eu E oo 
i _ 1] AT d z T — dE 
(so; Sa) E S (J Sa (w w tug 58a (v o 4) 5do 
w +B/2 
E Sa ) 5dw) 
= Sa (w-w +w )>Sa (w-w -w U) 
ESL Chiq 2 emd 2 
e 
HT oo 
d T T 
+e NN Sa (w+w 704) 58a (ua t0 4) 5,do 
-w +B/2 
É T T 
= f Sa (uto -w 4) 58a (wtw to.) 5do) 
-w -B/2 
E 
w +B/2 
pu T a 2 T 
- ( | sw-w *0.)zdu - f Sa (w-w -*u .)zdu) 
er veg e a 
-œ w -B/2 
C 
a SEE 
2 T S 2 T 
- ( | salwtw_ -w.)>dw - f Sa (uu -u .)zdu)] (4.9) 
cou Cr d 2 
«= 00 20 7B/2 


Equation 4.9 is certainly cumbersome, but both the 
number of variables and the complexity can be reduced by making 
the following change of variables on the positive frequency 


integrals: 
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l1) teta (w=w_)T/2 then (272) dx@= <din 


to, X = to 


2) when w : 


3) when w = w.tB/2, X = SEM Cm 


For the negative frequency integrals the following substitu- 


tions will be made: 


1) let x = (wtw )T/2, then (271) ax "da 
2) -whenu = 2o x = o 
3) when o = -w tB/2, xX = tBT/4. 


After performing these substitutions, it is seen that the 
first two integrals are equal in amplitude, but have opposite 


phase angles. As a result of this, (so Sq) becomes 


2 oo 
- 2 
(EE E RE Sa (x-o 7/2) dx 
BT/4 , 
E f Sa (x-04T/2) dx) 
-BT/4 
E BT/4 , 
el peo Zen c f Sa^ (x-u ,T/2) dx) 
—00 -BT/4 


00 


-2 cosugT( f Sa(x+w_¿T/2)Sa(x-w¿T/2) dx 


BT/4 


- f sa (x+w¿T/2) Sa (x-w T/2) dx (4.10) 
ZBD4 
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By similar analysis 


A'T 2 2 
(sis) = SIC f Sa(x-o,T/2)dx - Sa (X=) ,T/ 2) dx) 
ad 2 d -BT/4 d 
o 5 BT/4 , 
* ( f Sa(x-w,T/2)dx - Ji Sa cto ,7/2) dx) 
=% -BT/4 
E o | sa (x+w ,T/2) Sa (x-o 7/2) dx 
BT/4 
- f Sa (x+wsT/2) Sa (x-w¿T/2) dx) ] (4.11) 


-BT/4 


) 


The next term in the conditional mean is D 


which is given by Equation 2.18. As previously done, this 


equation is modified to become 


00 


l 


Lu = mae o EOS 

WwW +B/2 

= ON emu 

0) 7 dp J 
C 

-w +B/2 

- D S. (w)N. (-w) dw] (49902) 

-w —B/2 dp J 


Since H(w) is an ideal filter, in order to form the product 
QE ccu) and Sap e ? the spectrum of n; (t) must be found. 

The equation for 74 E) (see Equation 4.3) shows a truncated 
process similar to that found for the BPSK case. Thus N; (w) 


is determined to be 
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=3 (w-w, ) T/2 


O s E 
DR = YP 3/7 z[Sa(w w) 5€ 
-j (w+w] )T/2 
+ Sa (w+w] ) 5€ 
=3 (w-w,)T/2 uum 2 
a Sa (w=w,) 38 0 E Sa (u*u.) Fe 0 (a an 


Now, multiplying No (0) by S pt), where S, (o) is given by 


d dp 
Equation 4.9 yields eight terms. Four are on the positive 
side of the w= 0 axis and four are on the negative side of 


the w = 0 axis. For convenience just the positive frequency 


terms are shown here 


2 
= AT e T 
N LEN ETE "Pa; TlSa'to Wo) 5 
E: NEU I iu yl 
2 cosw ,TSa (w Wy) xSa (wtu,) > +Sa lw WF] (AMEN 


Expressing w, and Q4 aS u. tug as before, and substituting 


0 
Equation 4.14 into Equation 4.12 under the assumptions stated 
in the development of (si,Sq); yields an expression similar 

to Equation 4.9. The difference here is that the Sa (u-w tw 5) T/ 
2Sa (w-w wa) T/2 terms are multiplied by 2 cos wa T rather 

than ee By making the same change of variables as in 
Equation 4.10 and recognizing that Equation 4.14 covers only 


half of the frequency spectrum (which is symmetrical), 


(n',s.) becomes 
J] q 
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co 


uS =— 2 
(nj,SQ) = 2; AT P, /TII m Sa'(x-u ,T/2) dx 


sn M 
- f Sa“ (x-wqT/2)dx)( f sa” (x+w_T/2) 
-BT/4 E 


BA S 
- Í sá” (x+w T/2) dx) 
-BT/4 


co 


-2 cosuQT( fJ Sa(x+w¿T/2)Sa(x-w¿T/2)dx 


BT/4 
- f Sa (x+w¿T/2) Sa (x-w¿T/2) dx] (o) 
-BT/4 


| | » 1 2 2 
Finally, with (sj,s,), (nj,s,) and a (lsgll = eae ) 


determined, the conditional mean of G' can be expressed as 


2 
| : 7 EAT) ACD. 1 
E [G ER transmitted] = [AT/P. ;/T Ta] Is 


a o; Hu 
Mum sag 7/2) dx = f sa” (x+w¿T/2) dx) 
-o -BT/4 


a 5 BT/A , 
+ ( f Sa (x-04T/2)dx - f Sa (x-w4T/2)dx) 
-0 -BT/4 


co 


-2 cos o4T( f Sa (x-04T/2) Sa (x*»4T/2) dx 


= 00 


BT/4 
- f Sa (x-u,T/2)Sa(x*o ,T/2)dx)], i = 0,1 (4.16) 
-BT/4 
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2. Calculation of the Conditional Variance 
The conditional variance of G' in the frequency domain 
is given by Equation 2.27, which can be modified to take on 


the following form 


N 

VAR[G'|s. (t) transmitted] = E 
1 2T 

W +B/2 


[ f ER | cm - f E 
0 P w—B/2 


2 "hu 
ap? | doge i aom (4. L7) 


due to the fact that H(w) is assumed to be an ideal filter. 


Therefore in order to calculate the conditional variance, all 





that is needed is the determination of Sap tu) |2. Since 
2 . . } 
[S gp (0) | = o E and Sap «9! is given by Equation 4.8, 
the terms on the right hand side of the w = 0 axis are given by 
2m2 

p. a TAT 2 T2 T 
[Sag C2) | = à [Sa (o w tug) zsa lw m OF 
positive - I 
frequency -2cosw ,TSa (w-w_ tw 4) 58a (w-w 704)51 (4.18) 
terms s = 
where Wy = 0707 and w = UFU ae 


Since the limits of the integrals in Equation 4.17 are 
on the right hand side of the w = 0 axis, for large We and 
small B compared to Was the terms of [Sy (0) |^ on Ene” Pett eiame 
side of the w = 0 axis can be neglected. By substituting 
Equation 4.18 into 4.17 it is quickly recognized that the inte- 


grals are of the same form as Equations 4.9, 4.10 and 4.14. 
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So, instead of repeating the steps which have been demonstrated 


before, the conditional variance is shown here without further 


development 
AT B 
VAR[G' [s; (t) transmitted] = weite | Sa/(x-u 7/2) dx 
BT/4 oo 
- li Sa(x-w 1/2) dx) + ( n Sa{xtw4T/2) ax 
-BT/4 —o0 
BT/4 


Sa{xtw ,T/2) dx) 
-BT/4 


00 


-2cosw T( J Salx=w¿1/2)Sa (x+w¿T/2) dx 


BT/4 
- f Sa(x-w.T/2)Sa(x+w T/2)dx)], i = 0,1 (4.19) 
-BT/4 | 


3. Probability of Error Calculation 
With the conditional mean and variance of G' given by 
Equations 4.16 and 4.19, it is now possible to compute the 
probability of error performance of the BFSK coherent receiver 
with an ideal front-end filter. Assuming that the probability 
of the transmitter sending a "one" or a "zero" signal are 
equal (P (s, (£)) - P (sq (t)) = 1/2)), Equation 2.39 shows that 


the average probability of error P. is 


mv. 


o 2 2 
Pa = 3i eee aa y een adsl. (45205 
ERA / 2T = 0 /2T 
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where the threshold setting given by Equation 2.33 has been 
utilized. 

In order to compute Po’ all that is needed is the 
specific conditional mean divided by the square root of the 
conditional variance and using this ratio in the limit of the 
appropriate error function integral. Recognizing that 
A°T/2 = E, the average energy per bit, and making the appro- 


priate cancellations, ae is determined to be 


= /E/Ng (1 - /2Pn./E) YFilter Factor (4 21) 


| 
ls 
Q O 


and TA becomes 


re 


| 


E -VE/Nç (1 + (2P,/E) VFilter Factor (4.22) 


E 
a 


where the Filter Factor is defined as 


Filter Factor = E 
2T 
M Ed uw 
[( f Sa(x-w.T/2)dx - f sa(x-w.T/2) dx) 
d d 
-% -BT/ 4 
E m 
+ ( J Sa"(x-w5T/2)dx - f Sa(x=w¿T/2) dx) 


= -BT/4 


00 


- 2 cosw,T ( a Sa (x-wsT/2) Sa (x+w¿T/2) dx 


BT/4 
- | Sa (x-0¿7/2)8a (x+w¿T/2)dx)] (2:25) 
-BT/4 
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The probability of error performance can now be determined 
versus E/Ng, the signal to noise ratio; dB the jammer to 
signal ratio; and BT, the relationship of the filter bandwidth 
to the inverse of the signal bandwidth, by using Equations 
4.21 and 4.22 as the limits of integration in Equation 4.20. 
It is desirable to compare these results with those 
previously obtained for both the jamming problem and for 
the general FSK coherent receiver. In order to compare the 
results for the probability of error performance of the coher- 
ent FSK receiver with no filtering, it is necessary to analyze 
how those results were obtained. For the coherent receiver 
operating in just additive white Gaussian noise alone, the 


probability of error from Ref. l is 


p E J Elo au? /2 
= 


du (4.24) 


E 

Z 

© 
NO 
d 


This result is obtained when the energy per bit is A^T/2 
and the normalized cross correlation p is zero, where (see 


Ref. 1) 


f 
3 2 2 

E - el [yg (€) + yy (t)Jat (4.25) 

0 

È 
" fl f 
0 " g f Yolt)y, (t)dt (4.26) 
t 
0 


65 


and t t, are replaced by 0, T respectively. From these 


(ty 
expressions it is easy to show that wT must equal nm and 
WoT = kn (where n and k are integers) for E = A°T/2. For the 


FSK modulation p is 
p = Sa (w) +w) T + Sa (w, -wç)T (44.25 


For large (w,+tw))T the Sa(w,+w,))T is negligible and p is 

zero when (o4 709g) T is equal to mr (where m is an integer). 
When this is substituted into Equation 4.23, recalling that 
QT = (04709) T/2, it can be seen that when BT is zero the 
Filter Factor of Equation 4.23 is 1.0. For values of QT 

that are odd multiplies of 7/2 the Filter Factor is exactly 
equal to 1.0, but for even multiples it becomes 1 plus a small 
quantity that is negligible. So again, as in the PSK ideal 
front-end filter case, for BT equal to zero, Equations 4.53 
and 4.54 are in agreement with the results of Ref. 2 for the 
coherent FSK receiver operating against the optimum jammer. 
For Pnj/E equal to zero these results correspond to the results 


of Ref. 1 for the coherent FSK receiver operating in additive 


white Gaussian noise alone. 


B. FSK WITH A SECOND ORDER FRONT NDS PE te 

In order to evaluate the effect of a second order band- 
stop filter on the probability of error performance of a 
coherent FSK receiver, the filter specified by Equation 3.37 


will be utilized. The use of the equivalent baseband filter, 
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along with Parseval's theorem will not greatly aid in the 
determination of the conditional mean and variance of G' 

(the output of the receiver of Figure A.2). For convenience, 
the parameters of the coherent FSK receiver model are repeated 


here (see Equations 4.1 and 4.2) 


x m _ ed 
S, (t) = 2A sin 5 (wi wg) t sin 5 (w, +wp)t 
sg (t) AMOS Wot 
s} (t) = A cos wt UNE qu 


The optimum jammer derived in Ref. 2 is 

n.(t) = -yP ET =e Sn E diu rte O EL < T 
J DJT 2 NO) Dome © <0 — — 

and the filter of Equation 3.37 is 


- + JwB 
wg TY j 
where W 6 will be replaced by Da as in the FSK calculations 
involving the ideal filter case. 
l. Calculation of the Conditional Mean 
As before, Equation 2.7 will be used to determine the 
conditional mean of G'. Since the FSK signals are equal 


energy Signals, all that is needed from Equation 2.7 are the 
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terms Uu SM and SS ) with à = 0,1. These terms are given 


d 
by Equations 2.14 and 2.18. First (Sir Sg) will be calculated 
by evaluating the product of S; (-u) with Ee along with 


H (+w) and then substituting into Equation 208 which ae 


repeated here for convenience 


oO 


l : 
a) E 5, | J S; (u) Sgp VMC oldu i = 07I 


Investigating the particular case of (S9rSa) : the product 
ont Sg C) and Sap is given by Equation 4.8 where y _ Fw and 
Y_¿TW¿ were substituted for Oy and o, respectively. Thus, 


0 
substituting Equation 4.8 into Equation 2.18 yields 





(ss) = dl A^ T. 
Onc 27 4 
pias co 
d T T 

SES T saturu To RD 

E > T T 
e i Sa (ui o sa ee 
E í sao £9. )iH(-u)du - | suta -w 4) SH (-w) dw] (4.28) 


to further determine (SQ ( Sq) it is necessary to investigate 
both H(-w -2x/T) and H(w i -2x/T) . By substituting W 72X/T for 
w in Equation 3.37, expanding, and multiplying numerator and 


denominator by T/4u.. results in 
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H(».-2x/T) - E 


2 
(RE ao 


oc T 4 


(41:590) 


The filter transfer function was derived under the assumption 


that "M >> B so H (0 72x/T) becomes 


2 
EE 
o T 
H (w_-2x/T) = > (4.31) 
x JET 
E pt 4 
C 
and similarly 
x^ 
X t w T 
H(-w -2x/T) = Sore 15 2.) 
x SET 
ioo bem ay 
e 
For ease of notation let H4 (x) = H(w—2x/T) and H, (x) = H(-w -2x/T). 


Then separating H, (x) and H 


yields 


(x - 2)? 
Re{H, (x) } = 
x 2 BT 
T i 
2 
e dar 
Im(H,(x)] = | 
I x - 3? 4 (BE)? 
A oc T 4 


> (X) into real and imaginary parts 


(4.33) 


(4.34) 


(x "mn 
RetH, (x) ] = 5 (435) 
ee 
w T 4 
C 
2 
BT x 
pom T 
Im{H,(x)} = —s (4.36) 
x pr 2 
(x ‘oF po 


with Equations 4.33 through 4.36 it is now possible to expand 


Equation 4.29 into real and imaginary parts which yields 


2 E 
[NE “Gq [coswgT I f Sa(x-w¿T/2)Sa (xt ¿T/2) 


00 


[Re{H, (x) }+Re{H, (x) }dx 


00 


+ j E Sa (x-u ,T/2) 8a Go ,T/2) (Im(H, (x))eIm(B, (x) )) dx] 


00 


sinuo T Sa (x-w ,T/2) Sa (x+w,T/2) (Im{H, (x) }-Im{H, (x) }) dx 


+ 


+ J E Sa (X-u ,T/2) Sa (x *oQT/2) (Re(H, (x) J-Re(H, (x) J) dx] 


cO oo 


- [ 3 salta gT/2)Re(H (0) dx + f Sa‘(x-w 47/2) Re{H, (x) }ax 


00 


do) E Sa Gc 7/2) Im(H, (x) )Sa!(x-o 47/2) Im (B, (x) Jdx] ] (ARST) 
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With the real and imaginary parts of H4 (x) and H, (Xx) as given 


by Equations 4.33 through 4.36, it is easy to see that the 
imaginary terms integrate to zero because of the odd symmetry 
about x = 0 that the functions of Equation 4.37 exhibit. Thus 
1 
(SQ Sa) becomes 
-A^T 


(prsa) = al) Sa (xu 7/2) Re(H. (x) }ax 





t J Sa‘(xtw ,T/2) Re{H, (x) }ax 


- OO 


00 


- cosw 4T Sa (x-u ,/2) Sa («o 47/2) (Re(H, (x) }+Re{H, (x) J) dx 


=—O 


00 


- sinw¿T / Sa (x-w ,T/2) Sa (xtw5T/2) (Im{H, (x) }-ImtH, (x) }) dx] 


= © 


(4.38) 


and by similar analysis 


aT ue 
(SiS) = “Gol j Sa(xtogT/2) Re(H, (x) Jdx 


— CO 


+ f Sa x-u 47/2) Re(E., (x) )dx 


00 


- cosw¿T J Sa (x-w ,T/2) Sa (xt+w ,T/2) (RetH, (x) }+ Re{H, (x) }) dx 


00 


- sinw¿T _J Sa (x-w ,T/2) Sa (xtw ,T/2) (Im{H, (x) }-Im{H, (x) }dx] 
(2839) 
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The next term that must be computed is TS UR given 


by Equation 2.18 and repeated here for convenience 


oo 


em E = 
(ns,sg) = 3 ii Sgp (©) Nj (~w) H (~w) dw 


== OO 


As before, the calculation of this term starts with the pro- 


duct of ER, and S, (50). This product when fully expanded 


dp 


contains sixteen terms. Eight are products of Sa(x) functions 


centered on opposite sides of the u - 0 axis and, for large Ws 


Wo, are essentially equal to zero. By expressing Wy aS W Was 


and w, as w_+w substituting into Equation 2.18, and performing 


1l C 


the same change of variables as in the expansion of (Si rS)? 


it can be shown that mono becomes 


00 


(nissa) = PT f sá (to 1/2) (E (0 22x/T) +H (=u -2x/T)) dx 


= 00 


+ f Sa x-u,T/2) [H(w,-2x/T) +H(=w, -2x/T))dx 


— OO 


oo -JugT 
- 2 f Sa (x-w ,T/2) Sa (x+w T/2) [e H (w -2x/T) 


— OO 


juqQT 
+ e H (7o 7 2x/T) ) àx] (4.40) 


The integrals in Equation 4.40 are now of the same 


form as those in Equation 4.29. Therefore, the same change of 


12 


variables operation can be performed as done in the derivation 


of Equation 4.40. This results in 


Eus - Pal) Sa/(x-u ,T/2) (RetH, (x) }+RetH, (x) PNG 


oo 


+  f Sa(xtw 47/2) (Re{H, (x) }+Re{H, (x) }) dx 


=> 200 


oo 


- 2coso 4T J Sa (x-wsT/2) Sa (x+w T/2) (Re{H, (x) }+Re{H, (x) }) dx 


00 


T J Sa (x-0 ,T/2) Sa Go 47/2) (Im{H, (x) }-Im{H, (x) }dx] 


- 2sinw y 


(4.41) 


where Re(H, (x) 5, Re {H, (x)), Im{H, (x) } and Im{H. (x) } are defined 
by Equations 4.33 to 4.36. The conditional mean of G' under 
the assumption that Sp (t) or s (t) was transmitted can be 
expressed as the sum of Equation 4.38 or 4.39 with Equation 
4.41. 
2. Calculation of the Conditional Variance 
The conditional variance of G' is given by Equation 


2.27 and repeated here for convenience 


00 


N 
0 2 2 
VAR IG |s; (t) transmitted] = zm f [Say (0) | IH (o) | ^do 


— OO 


BD m ecognizing that HG) | ^ Powis,  Lomereal h(t), 


Equation 2.27 can be expressed as 


73 


00 


N 
| 0 
VAR[G'|s; (t) transmitted] = + Jf |s 


«== 00 


ap (©) |H (w)H (~w) du (4.42) 


In order to begin the calculation of the conditional variance, 
as for the ideal filter case, Sao (10) | must be determined by 
multiplying Ella by US This product contains sixteen 
terms as encountered in the product of Ni Aia] and E 
2 i 
Fortunately, ISag (0) | only differs from N) (=w) Sgp (v) by a 
constant. Therefore it is possible to perform a similar set 
of operations as done in the derivation of Equation 4.29 so 
as to arrive at 
NgA^T 
VAR[G'|s, transmitted] = 
i 4T 





[ E Sa/ix-u 1/2) [H (2x/T-0,)H (0 -2x/T) +H (-2x/T-0,)H (w,+2x/T)) dx 


+ 1 ii Salar T/2) (H(2x/T-u,.) H (o, -2x/T) *H (o, +2%/T) H (-w -2%/T) ) dx 


co Jw aT 
- 2 f Sa (x+w4T/2) Sa (x-w4T/2) (e 


H (2x/T-w JH (w 72x/T) 


-JugT 
+ e H (Zx/Tto )JH(-2x/T-w )) dx] i = 0,1 (4.43) 


where H (7o -2x/T) and H (o 72x/T) are given by Equations 4.31 
and 4.32. It can be demonstrated that H (2x/T-0.) and 


H (2x/Ttu.) are given by 
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x 
(x + TF! 
H (2x/Tłw ) = E o (4.44) 
(x 4 -—X- JBT 
WT 4 
and 
2 
(x - 37) 
H(2x/T-w.) = E O (4.45) 


By substituting Equations 4.30, 4.31, 4.44 and 4.45 into 
Equation 4.42 and expanding it into real and imaginary parts 


the conditional variance of G' becomes 


NgA^T 
' = 
VAR[G'|s; transmitted] = A 
2 2 
X x 2 
co ap poma. 
Í salx-w¿T/2) CARR RA + > 
=% x 2 BT, 2 x 2 BT. 2 
(x uos x (75-) (x top E (—) 
C E 
2 2 
X 2 X 2 

oo w gt o ar 

A 
=00 S BT, 2 X 2 BT, 2 
(x =p) ar 5 (X Borm T (—) 
C C 
2 2 
X 2 2 
oo w aT w aT oer es ta 7 
2200s) .T f Sa(x-——)S (xt——) [——————— ——9—— 4 —————À—9————————jdx 
d 2 D 2 
E e A BE iy) (BL) 
T 2 T 4 
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2 2 


AZ x 2 
co waT Wat = o T um 
-2cosw ,T i Sa (x-——) Sa (x+———) [-——_—_—_>—_—_——————_ 4 ——————————] dx 
d 2 2 2 2 
de AA E 

DT 2 O T 2 

E C 
i = 0,1 (4.47) 


3. Probability, of Errori Calcula cron 

The conditional mean of G' is given by the sum of 
Equation 4.41 with either Equation 4.38 or 4.39 depending 
whether Sg (t) Or s, (t) were transmitted, respectively. The 
conditional variance is Equation 4.47. With the conditional 
mean and variance determined it is now possible to compute the 
probability of error performance of the FSK coherent receiver 
with an ideal front-end filter. Assuming the probability oi 
the transmitter sending a "one" or a "zero" signal are equal 
(P (sy (t)) = P (sg (t)) = 1/2)), Equation 2.39 shows that the 


average probability of error P is 


g e 
d.e e fe ax) 


E (4.48) 


where the threshold setting given by Equation 2.33 has been 
utilized" 

In order to compute Poe all that is needed is the 
specific conditional mean divided by the square root of the 


conditional variance since this ratio must be used in the limit 
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of the appropriate error function.» However, the integrals of 
Equation 4.47 do not match those of the conditional mean, as 
found in every case investigated thus far. So before going 

on to compute the limits of the error functions of Equation 
4.48, the behavior of the conditional means and the variances 
will be investigated with respect to wT. Referring to Equa- 
tions 4.33 through 4.36 and 4.44 to 4.45, each has a term of 
the form x^/u 1 in both the numerator and denominator. 
Remembering that the filter was derived using an assumption 

of large Qr the effect of the x^/u, T terms can become negligi- 
ble when used in the integrals involving products of two Sa (x) 
functions, since these functions decrease rapidly with increas- 


ing x. Therefore for large oc T, a Filter Factor can be defined 


as 
oo waT w¿T x^ 

Filter Factor = f Sa (x+-5—) +Sa (x-—5— | (—5— 87-2) dx 

— OO x + (>=) 

2 
co Wat wT x^ 
- 2cosu,T f a ee Oa) ere oO (4.49) 
— OO xX + (5) 


Thus the conditional variance becomes 


A^TN 


var(G'|s, transmitted] = se ESE a CES (4.50) 





i = 0,1 


and the conditional mean becomes 
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i+] ACT 
E[G'|s, transmitted] = [ATVP, ¿77 =1 == 
x= x Filter Eactor i = 0,1 (4.519 


The Filter Factor was computed via numerical integration 

methods and compared with the complete expressions of Equation 

4.47 for oT as small as 50 and as large as 1000 with negligi- 

ble errors. An example of a large wT would be a communications 

system operating at 1.0 MHz with a bit length of 13.3 milli- 

seconds (75 baud). The wT term would have a value of 13,300 

which further supports dropping the x^/u, T term from the equations. 
Evaluating now mo/7v.. by dividing Equation 4.51 by 

the square root of Equation 4.50 where the Filter Factor is 


defined by Equation 4.49, mo/Y Yo becomes 


[A TAPER nj/T - ne T721 E TT a Factor 


+ ssa (4.52) 
YYg VA TN x À.- x Filter Factor 


Recognizing that A^T/2 = E, the average energy per bit, and 


substituting this in Equation 4.52 yields 


© 


= VE/Ng (1 - /2P,¿/8 /1/2-" xFilter Factor (4.5559 


E 
a“ | 


and similarly 


E 


= -= VE/No (1 + V2P_./E ¿/E) V1/27 xFilter Factor (4.54) 


E 
a 
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The probability of error performance can now be determined 
versus E/Ng, the signal to noise ratio; Pg/Ps the jammer 

to signal ratio; and BT (the relationship of the filter band- 
width to the inverse of the signal bandwidth), by using 
Equations 4.53 and 4.54 as the limits of integration in 
Equation 4.48. 

As with the ideal filter case, in order to compare the 
results for the probability of error performance of the 
coherent FSK receiver with no filtering, it is necessary to 
remember that those results were obtained using p = 0 and 
E = A^T/2, where p and E are defined by Equations 4.26 and 
4.25 respectively. Therefore, for large (w two) T, the 


Sa(w,tw,)T term is negligible and p is zero when (0,-07)T is 


0 
equal tom . When this is substituted in Equation 4.23, 


recalling that w T = (04709) T/2, it can be seen that when BT 


d 


is zero the Filter Factor is 27. For values of WT that are 

odd multiples of 1/2 the Filter Factor is exactly equal to 

27, but for even multiples it is equal to 2 plus a quantity 
which is negligible. So again, as in the FSK ideal front-end 
filter case, for BT equal to zero, Equations 4.53 and 4.54 are 
in agreement with the results of Ref. 2 for the coherent FSK 
receiver operating against the optimum jammer. For EE 

equal to zero, these results correspond to the results of Ref. 

l for the coherent FSK receiver operating in additive white 
Gaussian noise alone. Again, as with the coherent BPSK receiver, 


as the BT approaches zero, the behavior of the probability of 


error performance of the coherent FSK receiver calculated using 
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a second order front-end filter is similar to the performance 
calculated using the ideal front-end filter. Later the per- 
formance comparison will be performed throughout the range of 


pertinent values for the variable parameters. 
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y" RESULTSORAND DENSCUSSTON 


A. GENERAL 

For each of the four receiver configurations (PSK and FSK 
receivers with both ideal and second order front-end filters), 
the results for the probability of error performance have been 
calculated and plotted as functions of D (jammer to signal 
ratio), E/Ng (signal to noise ratio) and BT (the filter band- 
width, bit duration product). The meaning of the jammer to 
signal and signal to noise ratios are self explanatory. 
However, the term BT deserves some attention. Since T is the 
time duration of sy (t) and So (tt), the spectra of each have 
been shown to be functions of the form Sa(uT/2) where Sa(x) 
has been defined as sin(x)/x. The Sa(wT/2) function has nulls 
when the value of wT/2 is an integer multiple of m or when 
IT). 

From the Sa(uT/2) spectrum it can be shown that the value 
of B = 27/T represents approximately the 3 dB (or half power) 
bandwidth of the signal (actually 27/T is the 4 dB bandwidth). 
The term BT can thus be thought of as the ratio of the filter 
bandwidth to the signal bandwidth for intuitive purposes. It 
is still convenient to refer in further discussions to the 
location of the nulls of the Sa(uT/2) function. 

The analysis of the calculated probability of error per- 
formance for each receiver configuration uses three types of 


plots. The first is the familiar plot of probability of 
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error versus signal to noise ratio, with constant values of 
BT and CEE The values of “q used are 0, .1, .5, 1, and 
l0. These values correspond to no jammer present, and jammer 
to noise ratios of -10 dB, -3 dB, 0 dB, and 10 dB. The par- 
ticular values of BT depend on the modulation structure 
analyzed, and will be discussed later. The second type of 
plot will allow the analysis of the probability of error 
performance to be viewed as a direct function of BT. That 
is, for a constant value of signal to noise ratio, the proba- 
bility of error is plotted versus BT. This plot consists of 
a family of curves corresponding to values of jammer to signal 
ratio (J/S) of 0.1 (=10 GB), 0:25 (=6 GB). 0-5) (CS amam 

1 (0 dB). Finally the third type of plot used in this analy- 
Sis is again a probability of error versus signal to noise 
ratio plot. However this plot will consist of a family of 
curves for values of BT appropriate to the case analyzed. In 
these plots the value of J/S used is zero. This corresponds 
to the case where the jammer waveform is considered not to be 
time truncated. 

The case J/S = 0 represents an alternate to the assumption 
that Be) is a time truncated process. When n; (t) is a pure 
cosine function, SE is a delta function, which can be 
completely removed by the front-end filter. The resulting 
plots then describe the effect of the signal being distorted 
by the front-end filter before it is correlated with Sq(t). 


The probability of error curves for this analysis is shown by 
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the set of plots described in the preceding paragraph. This 
performance analysis can also be viewed in each of the other 


plots with the curve J/S = 0. 


B. PHASE SHIFT KEY RECEIVER ANALYSIS 

l. General 

The probability of error performance of the BPSK receiver 

with both the ideal and second order front-end filters is 
anal” zea and compared to that of the unfiltered BPSK receiver 
in this section. The results of Equations 3.32 and 3.33 are 
substituted into Equation 3.29, yielding the average proba- 
bility of error for the BPSK receiver with an ideal front-end 
miter. Substituting équations 3.69 and 3./0 into Equation 
3.66 yields the average probability of error for the BPSK 
receiver with a second order front-end filter. In each case 
the plots described in the beginning section of this chapter 
are used to analyze the performance. For PSK modulated sig- 
nals, the discrete values of BT are chosen with respect to 
the PSK modulation frequency spectrum. Recalling that the 
spectrum (Equation 3.22) is of Sa(wT/2) form, the first null 
of this spectrum occurs for w = 27/T. Thus the null to null 
bandwidth is 4m/T. Therefore, the value of BT = 47 will be 
used as a practical limit in this analysis. The other dis- 
crete values of BT are: BT = 0 (no filter), BT = 0.4 (a 
small notch), BT = 3.16 (25% ot the null to null bandwidth), 
BT = 6.28 (503% of Ene null to null bandwidth), and BT = 12.56 


(the full null to null bandwidth). 
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2. PSK Receiver With An Ideal Front-End Filter 

Figures A.3 through A.9 result when the probability of 
error for the BPSK receiver with a front-end filter is calcu- 
lated using the ideal filter model. These plots show how the 
probability of error performance is affected by the signal to 
noise ratio, jammer to signal ratio and BT. Figure A.3 repre- 
sents the probability of error performance of the BPSK receiver 
with no front-end filter. This corresponds to the results 
derived and plotted in Ref. 2. In fact the J/S = 0 curve in 
this plot is the well-known result for the coherent BPSK re- 
ceiver operating in additive white Gaussian noise alone. This 
curve will serve as a performance reference in the following 
discussions. Figures A.3 through A.8 allow the probability of 
error performance to be graphically analyzed for the cases 
where the jammer is either a truncated precess of a continuous 
time function as discussed in the first section of this chapter. 

For the truncated jammer, by following the J/S = 0.1 
and 0.5 curves from Figure A.3 to Figure A.8, the curves shift 
to the right as the value of BT is increased. This means that 
for a given probability of error, the receiver operating with 
the larger value of BT requires more signal to noise ratio 
to achieve a given level of performance. Thus for the optimum 
jammer (truncated by the receiver process) the front-end filter 
actually further degrades the probability of error performance 
for the BPSK receiver. Specifically, (refer to the J/S = 0 
curve of Figure A.3) for BT = 3.16 there is 5 dB of degradation, 


for BT = 6.28 there is 6 dB of degradation and for BT = 12.56 
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there is 9 dB of degradation. The rapid increase of the 
probability of error as BT increases is seen in Figure A.8. 
For a fixed value of E/N = 100 (20 dB), the probability of 
error rapidly increases as BT approaches 47. As BT increases 
further, the slope of the probability of error curves become 
less positive as less energy is incrementally removed from the 
signal by the filter. 

Assuming now that the jammer is not truncated by the 
receiver process, the analysis of the probability of error 
curves is performed by comparing the J/S = 0 curves of Figures 
A.4 to A.8 with that of Figure A.3. These are plotted together 
in Figure A.9. From Figure A.3, for J/S equal to -10 dB, 

5 dB of additional signal to noise ratio is required to achieve 
the same probability of error performance as with no jammer 
present. For J/S equal to -3 dB, 10 dB of additional signal 
to noise ratio is required, and for J/S greater than one, no 
increase in signal to noise ratio results in improved perfor- 
mance as discussed in Ref. 2. Assuming that the filter com- 
pletely removes the jammer, Figure A.9 shows the performance 
penalty associated with the use of the front-end filter. By 
comparing Figure A.9 with Figure A.3, it is seen that for 
large values of J/S the receiver probability of error perfor- 
mance is improved by the front-end filter. The specific level 
of improvement is determined by comparing the curve in Figure 
A.3 corresponding to a particular J/S value with the curve in 
Figure A.9 for a fixed value of BT. For the ideal jammer in 


which the frequency spectrum is a delta function, an infinitely 


85 


narrow bandwidth filter can be used to counter this jammer. 
In a more practical sense, the filter bandwidth should be made 
as narrow as possible, and yet completely null the jammer over 
its entire spectrum. The sensitivity of the probability of 
error performance to the value of BT used for this case is 
shown in Figure A.8 with the J/S = 0 curve. Again, there is 
a rapid increase of probability of error as BT increases to 
the value of 12.56. This further shows that BT should be kept 
as small as practically possible. 

3. PSK With A Second Order Front-End Filter 

Figures A.10 through A.16 result when the probability 

of error for a coherent BPSK receiver with a second order 
front-end filter is plotted versus E/N, J/S, and BT. The 
probability of error in this case is obtained by substituting 
Equations 3.69 and 3.70 nto Equation 3.66. Figure A.10 
shows the probability of error performance with BT = 0 (actually 
BT = 0.0001 for computational purposes). For values of BT 
less than 3.16 (approximately equal to m), the probability of 
error curves compare closely with the curves plotted for the 
BPSK receiver with the ideal filter. This is illustrated by 
comparing Figure A.ll with Figure A.4, and Figure A.12 with 
Figure A.5. For values of BT greater than 27 the comparison 
of Figures A.6 and A.7, with Figures A.13 and A.14, show that 
for BT = 6.28 the ideal filter calculations differ from the 
second order filter calculations by 1 dB. For BT = 12.56 


there is a 2 dB difference in the curves calculated using the 
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ideal versus the second order filter. The comparison of 

Figures A.8 and A.15 show that as BT increases, the calculations 
using the ideal filter are less accurate compared to the second 
order filter calculations. Again, for the truncated jammer 
case, no value of BT results in improved probability of error 
performance in comparison to the unfiltered receiver. 

Now, assuming that the receiver process does not trun- 
cate the jammer, Figure A.16 shows the probability of error 
performance achieved for BT = 0, 3.16, 6.28 and 12.56. Com- 
paring Figure A.16 to Figure A.9 shows that there is similarity 
in the probability of error performance calculated using the 
ideal filter model with that calculated using the second order 
filter model, for values of BT less than 3.16. Therefore for 
values of BT less than m (or B less than 25% of the signal null 
to null bandwidth), the probability of error can be determined 
from the results derived for the ideal filter model being used 


as a front-end filter. 


EN CONERENTI FREQUENCY SAIPFT KEY RECEIVER ANALYSIS 

The probability of error performance of the FSK receiver 
with the two filters, the ideal and second order front-end 
filters, is analyzed and compared to the unfiltered coherent 
FSK receiver performance in this section. The results of 
Equations 4.53 and 4.54 are substituted into Equation 4.48, 
yielding the average probability of error of the FSK receiver 
with an ideal front-end filter. Substituting Equations 2.28 


and 2.29 into Equation 2.23 yields the average probability of 


87 


error for the FSK receiver with a second order front-end filter. 
In each case the plots described in the beginning section of 
this chapter are used to analyze performance. For FSK modu- 
lated signals, not only must the parameter BT be chosen, but 
also the frequency separation of the FSK signals must be set. 
From the general results on performance for the coherent re- 
ceiver, the probability of error is given by (from Ref. 1 


assuming no filters and no jamming present) 


e (1-p) 





P, = erfel (Sey 
/Ng£ (179) 
with 
n 
erfc(x) = f e / du 
x Y 2T 


where c and p have been defined by Equations 4.24 and 4.26, 
respectively. The results of Ref. 2 given by Equation 4.5, 
with the Filter Factor equal to 2 , use a value of p of zero. 
Therefore the comparisons made here will be for values of 

Q 170g such that p is zero. The term used to describe this 
separation and used in the derivation of the probability of 
error results is WaT, defined in Equation 4.7 as 4T = (04709) T/2. 
The minimum value of QaT is chosen to be 7/2 because from 
Equation 4.27 it is seen that the separation (04 70g = 71) “rs 


the minimum separation required for the normalized correlation 


coefficient of Sp (t) and S} (t) to be zero. The other value of 


88 


waT chosen is 1, where again with this separation the normalized 
correlation coefficient is zero. 
In Section IV.A.1 the strategy of placing the filter at 


the midpoint of wi and Wa is described. From Equation 4.13, 


0 
it is seen that the optimum jammer for the coherent FSK receiver 
consists of Sal T/2) functions centered at 01 and Wo - There- 
fore in order to filter out the major portions of the jammer, 
the bandwidth of the filter must be at least equal to the fre- 
quency separation of the FSK signals. This also holds true 
for the case where the jammer is considered not to be truncated 
by the receiver process. Given these conditions, for QT equal 
to 1/2 the values of BT used will be Zero, corresponding to 
no filter, mt corresponding to the difference of the signal 
frequencies, and finally m + 0.171 (= 3.46) in order to take 
into account the non-truncated jammer having non-zero spectral 
width. Similarly for WaT — m, the values of BT chosen are 
BT = 0, BT = 2 , and BT - 2m * 0.1:2m- (» 6.91). 
l. FSK Receiver With An Ideal Front-End Filter 

Figures A.17 to A.24 result when the probability of 
error for the coherent FSK receiver with a single ideal front- 
end filter is plotted as functions of signal to noise ratio, 
jammer to signal ratio, frequency separation, and BT. Figures 
A.17 and A.20 present the probability of error for the coherent 
FSK receiver operating with no front-end filter. Also the 
J/S = 0 curve in these figures correspond to the probability 


of error performance for the FSK coherent receiver (with 
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orthogonal signals) operating in white Guassian noise alone. 
It is not surprising that Figures A.17 and A.20 are exact 
duplicates, because in each case the normalized correlation 
coefficient of s,(t) and sg (t) is zero, as described in 
Section V.C. As with the BPSK analysis, the J/S = 0 curves 
of Figures A.17 and A.20 will serve as a performance reference. 

Assuming the truncated version of the jammer, first 
for W 4T = 1/2, Figures A.17 to A.19 show that as BT increases 
from O to 3.46, the probability of error curves shift to the 
right. This indicates that for increasing filter bandwidths, 
more signal to noise ratio is required to achieve the same 
probability of error performance as with no filter. From 
Figure A.17, for a value of BT equal to 1, approximately 2 
dB of additional signal to noise ratio is required to achieve 
the same level of performance as in Figure A.17. This applies 
for all values of J/S. From Figure A.l9, 2.8 dB of additional 
signal to noise ratio is required to achieve the performance 
level presented in Figure A.17. The increase in probability 
of error with increasing BT for all values of J/S is shown in 
Figure A.23. 

With wat equal to v, more filter bandwidth is required 
in order to attempt to notch out the jammer. As a result, a 
greater amount of the signal energy is removed by the filter 
than is the case for the receiver operating with wgT oe 2 
Therefore, the performance of the FSK receiver with a front- 


end filter and w,T set to the value of m, will be worse. 


d 


Figures A.20 to A.22 show that for BT increasing from zero to 
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6.91, more signal to noise ratio is required to achieve the 
performance level of the receiver operating with an unfiltered 
jammer. Specifically, 5.6 dB of additional E/N is required 
for BT = 21, and 6.6 dB of additional E/N is required for 
BT = 6.91. Figure A.24 shows the sensitivity of the probability 
of error performance for increasing values of BT. Focusing 
on Figures A.17 to A.24, in either case (wgT = m/2 or 1.57), 
analysis using the truncated jammer assumption shows that no 
value of BT causes a decrease in the probability of error for 
the coherent FSK receiver. Only degradation is observed. 

The analysis of the coherent FSK receiver with the 
deal front-end filter, assuming the Tion—truncated jammer, is 
performed by comparing the J/S = 0 curves on Figures A.17, 
Ae and A.l9 for w QT = 1/2. For wT = "7, the J/S = 0 curves 
on Figures A.20, A.20 and A.22 are compared. As before, the 
unfiltered J/S = 0 curves (Figures A.17 and A.20) serve as a 
reference. As in the BPSK case, this analysis shows that for 
either o 4T = 1/2 or m, the J/S = 0 curves still shift to the 
right as BT is increased indicating that more signal to noise 
ratio is required in order to achieve the same performance as 
for the unfiltered case. From the curves for an unfiltered 
jammer (Figures A.17 and A.20) for values of J/S greater than 
0.5, Lhe receiver is essentially inoperable (P, = 0.5). However, 
the J/S = 0 curve of Figure A.18 shows that an increase of 
E/N, of D OmebEcaneoning the performance of the filtered FSK 
receiver to the same value of p as for the unfiltered receiver. 


Figure A.18 corresponds to o SH cp o—cuPIigure A.lI9 


at 
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shows that for BT = 3.46, 2.8 dB of additional E/Ng is required 
to achieve the unfiltered, unjammed performance. Figures A.18 
and A.19 show that for values of J/S greater than .5 (for the 
coherent receiver) the use of the ideal filter does show an 
improvement on the performance of the FSK receiver when the 
filter is used to notch out the cw jammer. The sensitivity 

of the probability of error to the value of BT used for the 


case of w = 1/2, is shown by the J/S = 0 curve of Figure 


at 
A.23. As in the BPSK receiver case, this shows that the smallest 
possible filter bandwidth should be utilized. For Wgl = T, 
similar results are shown in Figures A.21 and A.22. With 

BT = 27, 5.6 dB of additional E/No is required to achieve the 
performance level of the unfiltered, unjammed receiver, and 

for BT =-6.91; 6.61dB of additional E/No 1s required. Figure 
A.24 shows the sensitivity of the probability of error for 

this case. An important fact to remember is that when using 
Figures A.23 and A.24, the ideal filter does not alter the 


jammer until BT is at least 2xo4T. 


2. The FSK Coherent Receiver With A Second Order 
Front-End Filter 


Figures A.25 to A.32 result when Equations 4.53 and 
4.54 are substituted into Equation 4.48 and are plotted versus 
BT, J/S and E/Ng. Figures A.25 and A.28 are the probability 
of error curves for the unfiltered FSK coherent receiver 


Operating with o,T = 1/2 and m, respectively, and with an 


d 
unfiltered jammer. As described In the analysis ones 


A.17 and A.20, Figures A.25 and A.28 are exact duplicates 


OZ 





because p in each case is zero. In general, the curves plotted 
using the second order filter model match within 1 dB to those 
calculated using the ideal filter model. Specifically, Figures 
A.26 and A.27 (wT = 1/2, BT = T and 3.46 respectively) 
compare within 1 dB to the results of Figures A.18 and A.19 
which are calculated using the ideal filter model. Similarly, 
[ET o QT = T, Figures A.29 and A.30 compare within 1 dB to the 
results presented in Figures A.21 and A.22 respectively. The 
sensitivity analysis from Figures A.31 and A.32 (o 4T = 1/2 

and m) still shows the rapid increase in probability of error 
with increasing BT. However, compared to similar curves calcu- 
lated using the ideal filter, Figures A.31 and A.32 have 
smaller (positive) slope. In all, assuming the truncated 
jammer, Figures A.25 to A.32 still show that as the front-end 
filter bandwidth increases, the probability of error increases 
also. Therefore, as with the ideal filter case, the second 
order front-end filter only further degrades the performance of 
the coherent FSK receiver operating against the truncated 
optimum jammer. 

Assuming the non-truncated version of the optimum 
jammer, it is easily seen that a second order front-end filter 
such as the one given by Equation 3.37 would not be completely 
effective against that jammer. By expanding Equation 4.3, 
it is seen that the jammer spectrum would consist of delta 
functions centered at w, andw,. The filter described in 


i 0 


Section III.B has a null only at w = two’ Therefore the 
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analysis cannot be performed by simply letting J/S = 0 in 
this case. From the analysis completed for the ideal filter 
case, it is seen that a higher order filter, or a filter with 


nulls at Wy and w should be used to defeat the non-truncated 


0’ 


optimum jammer. 
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VI. CONCLUSIONS 


A. GENERAL CONCLUSIONS 

In this thesis the effectiveness of front-end filtering 
techniques has. been investigated as an electronic counter- 
countermeasure (ECCM) in order for the coherent digital 
receiver to be able to operate effectively against the optimum 
jammer derived in Ref. 2. The specific coherent digital re- 
ceivers analyzed were the coherent BPSK receiver and the 
coherent BFSK receiver. The filters placed in the front-end 
Of the receivers, as shown in Figure A.2, were an ideal filter 
and a second order filter of the form given by Equation 3.37. 

Because the analysis of Ref. 2 defines the optimum jammer 
in the time interval (0,T), the analysis was performed using 
this interval in the calculation of the Fourier transforms 
of the jammer models. As a result, Sections V.B.2 and V.B.3 
as well as Sections V.C.2 and V.C.3 show that none of the 
front-end filters resulted in any decrease in the probability 
of error in comparison to receivers operating without front- 
end filters for the class of receivers analyzed. This occurs 
because the resultant jammer consists of a signal that is of 
the same form as sq), the correlation signal. As a result, 
the filter affects equal portions of both the signal spectrum 
as well as the jammer spectrum. Therefore under this analysis, 


the front-end filters further added to the probability of error 
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caused by the jammer and the additive noise for both the BPSK 
and BFSK coherent receivers. 

Taking into account the physical aspect of this problem, 
an ambiguity is seen. Both S(t) and Dos are defined over 
the interval (0,T). However over the next interval (T,2T) 
these signals remain unchanged. This repeats itself for the 
next interval of length T, and so on. In other words, both 
sq) and a are cw signals. Assuming such a form for 
the optimum jammer, the analysis was carried out by substi- 
tuting 0 for the value of J/S, for all but the FSK receiver 
with the second order front-end filter. Setting J/S = 0 pro- 
vides the same result as setting the term EPI equal to 
zero. In the case of the FSK coherent receiver, the second 
order filter does not have nulls at the frequencies of the 
Optimum jammer. This analysis cannot be performed simply by 
setting J/S = 0. 

When the PSK coherent receiver was analyzed with both 
filter models using the cw form of the optimum jammer, the 
filter was observed to introduce improvements to the proba- 
bility of error performance provided that the receiver was 
jammed with a value of J/S of at least 0.5. This is quite 
apparent for J/S equal to 1.0, as the unfiltered PSK coherent 
receiver becomes essentially inoperative (P = 0.5). When 
the filter 1s introduced in this case, the jammer term 
eS ae is set to zero. Therefore the conditional mean of 


G', the output of thessecelver isa function of (si,s) only. 
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When this happens, any increase in the probability of error 
compared to that of the unfiltered, unjammed receiver is 
caused only by the filter distorting the signal. The use of 
the filter is not without cost. As shown in Figures A.9 and 
A.16, additional signal to noise ratio is required to achieve 
the same probability of error as that of unfiltered unjammed 
receiver, 

The analysis of the FSK coherent receiver with the ideal 
front-end filter, under the assumption of the cw jammer shows 
similar results. However, because a single filter, centered 
at the midpoint frequency between the FSK signals is used, 
more signal to noise ratio is required to achieve the same 
probability of error performance as that of unfiltered, un- 
jammed receiver. This results because the filter bandwidth 
must be at least equal to the frequency separation of the FSK 
signals. Due to this observation, the smallest practical fre- 
quency separation should be used. Analysis of the second order 
front-end filter shows that this type of filter centered midway 
between the FSK signals was not as useful as a filter with 


nulls at the jammer frequencies. 


Be SPECIFIC CONCLUSIONS 
Some specific conclusions follow from the analysis per- 
formed in this thesis. These are as follows: 
1) For the PSK coherent receiver with a front-end filter, 


the probability of error performance calculated using 


9 


the ideal front-end filter matches the performance 
calculated using the second order filter within 1 dB 

for values of BT of less than m. Therefore, the analysis 
can be performed using the ideal filter for small values 
of BT. So, though this will require computer integra- 
tion methods, time will be saved by not analyzing complex 
filter models. 

Under the assumption of a time truncated jammer, the 
probability of error for the FSK coherent receiver 
calculated using the ideal filter compares with the 
results calculated using the second order filter with 

BE = 2wgT. Therefore, the ideal filter can be used for 
these calculations. 

The second order front-end filter centered between 

the FSK signals does not effectively remove the jammer 
(assuming the non-truncated jammer) and therefore is 

not as useful as a filter with nulls at the jammer 
frequencies. 

It appears that part of the performance degradation 
introduced by filtering is due to the fact that the 
filter distorts the signal components S(t) and Sp tt), 
so that the correlator receiver (or equivalently the 
matched filter receiver) is not matched to the incoming 
signals. This degradation could be overcome by using a 
matched filter receiver that is matched to s, (t) and 

so (t). While this has not been inves Ear amas 
thesis, it appears that such investigations would be of 


great interest and worthwhile. 
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Figure A.3. 


Coherent PSK Receiverwith an Ideal 


BT = 0.0 
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